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ABSTRACT 



Context, More than one hundred GeV pulsars have been detected up to now by the LAT telescope on the Fermi gamma-ray ob- 
servatory, showing peak energies around a few GeV. Current modelling proposes that the high energy emission comes from outer 
magnetospheric gaps, however radiation from the equatorial current sheet which separates the two magnetic hemispheres outside the 
light cylinder has also been investigated. 

Aims. In this paper we discuss the region right outside the light cylinder, or "near wind" zone. We investigate the possibility that 
synchrotron radiation emitted by thermal populations in the equatorial current sheet of the pulsar wind in this region can explain the 
lightcurves and spectra observed by Fermi/LAT. 

Methods. We use analytical estimates as well as detailed numerical computation to calculate the y-ray luminosities, lightcurves and 
spectra of y-ray pulsars. 

Results. Many of the characteristics of the y-ray pulsars observed by Fermi/LAT can be reproduced by our model, most notably 
the position of these objects in the P — P diagram, and the range of y-ray luminosities. A testable result is a sub-exponential cutoff 
with an index b = 0.35. We also predict the existence of a population of pulsars with cutoff energies in the MeV range. These have 
systematically lower spindown luminosities than the Fermi/LAT detected pulsars. 

Conclusions. It is possible for relativistic populations of electrons and positrons in the current sheet of a pulsar's wind right outside 
the light cylinder to emit synchrotron radiation that peaks in the sub-GeV to GeV regime, with y-ray efficiencies similar to those 
observed for the Fermi/LAT pulsars. 



Key words. (Stars:) pulsars: general - Radiation mechanisms: non-thermal - 
processes 



Stars: winds, outflows - Gamma rays: stars - Relativistic 



1. Introduction 



Since the launch of the Fermi space telescope, the sample of 
gamma ray pulsars has grown to include more than one hun- 
dred objects. All these pulsars, whether young or millisecond, 
single or in binaries, show spectra consistent with power laws 
with exponential cutoffs, where the cutoff energy lies in the 
range 1-10 GeV dAbdoetalJ l2010h . A prominent exception 
is the Crab pulsar, which has been detected by ground-based 
Cerenkov arrays in the TeV regime, and the s pectrum of which 
can be well fitted by a broken power law ( A iu et"aD l2008t 



IMcCann & for the VERITAS Collaborationl201 ll) . These obser- 



vations point to the outer gap/slot gap models as the most prob- 
able explanation for the pulsed gamma-ray emission, based on 



(Gonthieret al. 20101 


Decesaretal. 201 lb Pierbattista et al. 


2011 


: Venter et alJ 2009 


: Watters & Romani 201 ll Venter et al. 


2012 


), although for millisecond pulsars one has to evoke non- 



dipolar field geometries or displaced polar caps, in ord er to reach 
the required energies dHarding & Muslimovl 1201 ll) . In those 
models, the emission originates within the light cylinder, defined 
by the cylindrical radius r^c = c/io from the pulsar's rotational 
axis, where a corotating partic le would reach the speed of light. 
However it was pointed out bv lBai & Spitkovskvl(l2010l) that the 
emission region must be extended slightly outside the light cylin- 
der in order to reproduce the doubly-peaked light curves using 



the magnetic field configuration from self-consistent, force-free 
simulations of pulsar magnetospheres. 

The idea that high energy pulsations might come from 
current sheets ne ar or outside the light cylinder is not new. 
iLvubarskiil (1 19961) predicted that particles accelerated though re- 
connection close to the light cylinder might emit gamma-rays 
through the synchrotron process. Lyubarsky's emission site is 
the point where the warped equatorial current sheet which sep- 
arates regions of opposite magnetic field polarity in the pulsar's 
wind meets the curren t flowing to/from t he pulsar's polar caps, 
the so-called Y-point (ISpifkovskvll2006l) . Later it was pointed 
out that pulsed emission naturally arises from the periodicity of 
the pulsar wind, which is modulated by the star's rotation, in 
combination with the large bulk Lorentz factors of the outflow 
et al.ll2002l) . This model has been successfully applied to 
calc ulate the polarizatio n of optical emission from the Crab pul- 
sar dPetri & Kirk! 120051) and it has also shown promise in ex- 
plaining the gamma-ray light curves observed by Fermi/LAT 
(lPetrill2011l) . In these models the emission is attributed to the 
inverse Compton process and starts in the wind region far from 
the light cylinder, r » ric- In the same context it was proposed 
that the gamma-ray radiation of the Fermi/LAT band could be 
synchrotron radiation from power-law electrons in the pulsar 
wind ' s equatoria l sheet far from the light cylinder dKirk et al.l 
120021: lPetrill2012h . 

There is, however, no obvious reason why emission should 
be truncated outside the light cylinder (for magnetospheric gap 



1 



I. Arka & G. Dubus: Pulsed y-rays from thermal particles in pulsar wind current sheets 



models) or start at a minimum radius in the far wind zone (for 
the striped wind model). The emission region should evolve 
continuously from the outer magnetospheric gaps through the 
light cylinder and into the equatorial current sheet, possibly af- 
fected by the reconnection proces s that will inevitably occur in 
the sheet dBai & Spifkovskvll2010h . The region beyond the light 
cylinder but still in the "near wind" zone is the area that we 
are trying to explore. The purpose of this paper is to demon- 
strate how emission in the 100 MeV-100 GeV range explored 
by Fermi/LAT can naturally arise from thermal populations of 
particles in the striped wind, without going into the details of 
reconnection physics, but just using a few simple assumptions 
about the local description of the current sheet. In section 2 we 
describe the model we use, including our assumptions, and give 
some analytical estimates for the emitted radiation. In section 
3 we present examples of emission maps and phase-averaged 
spectra which arise from our model and explore the parameter 
space in order to come to more general conclusions about the 
pulsar population observed by ferm//LAT. Finally we discuss 
our results and the potential of a more detailed description of 
the current sheet outside the light cylinder as a source of pulsed 
energetic radiation. 

2. Equatorial current sheet 

2.1. Description of the particle population 

In order to describe the wind outside the light cylinder, we will 
use the "slow-rotator" solution to the magnetohydrodynamics 
equation s governing the ph ysics of the pulsar wind, which was 
found bv lBogovalovldl999l) . This solution refers to an electron- 
positron wind launched by a rotating neutron star, the magnetic 
axis of which is at an angle x to its rotational axis. In this de- 
scription the wind is purely radial and super-fast magnetosonic 
at launch, and the field can be described by a radial and an az- 
imuthal component of respective magnitude: 



B r 

Bu, 



Blc 
R 2 
Blc 
BR 



sin# 



(1) 
(2) 



where R = r/r^c is the spherical radius normalized to the light 
cylinder, § is the polar angle, B — (1 - 1/T 2 ) 1 ^ 2 is the bulk speed 
of the wind, normalized to the speed of light and Blc is a fidu- 
cial magnetic field magnitude at the light cylinder. The adequacy 
of this solution in describing approximately the structure of the 
wind for radii as small as the light cylinder has bee n confirmed 
by force-free simulations of pu lsar magnetospheres dSpifkovskvl 
1200(1 iBai & Spitkovskvl20Ta) . 

A prominent feature of Bogovalov's solution is the warped 
equatorial current sheet, separating the two magnetic hemi- 
spheres, across which the field changes sign. In the mathemati- 
cal solution the sheet is just a discontinuity, however in reality 
it should have a finite thickness. Such a current sheet is popu- 
lated by hot particles, the pressure of which balances the mag- 
netic pressure of the cold, strongly magnetized plasma outside 
the sheet. The sheet oscillates in space and time with a wave- 
length of A — 2nB/oi, and with the pulsar's frequency 10 - 2n/P, 
where P is the pulsar's period measured in seconds. 

We will assume that, in the wind frame (the frame that prop- 
agates radially outwards with a Lorentz factor equal to the bulk 
Lorentz factor of the wind, F), the shee t can be loc ally described 
by the relativistic Harris equilibrium (Hohl ll966l) . For this de- 
scription to be valid, the segment of the sheet under considera- 



tion should be approximately flat. For a relativistic outflow, the 
hydrodynamically causally connected region of the flow has an 
opening ange of roughly ~ 1/F, centered on the direction of 
motion of the outflow. Since for a pulsar wind F » 1, the hydro- 
dynamically connected sheet segment can be considered flat and 
the local Harris sheet description should be a good approxima- 
tion. 

We will denote quantities measured in the wind frame with 
a prime. If we denote the sheet normal direction in the wind 
frame with a capital X', and the sheet midplane is at X' = 0, then 
locally the field inside and outside the segment can be described 
by a tangent hyperbolic profile: 



B'(X') = ±B' tanh 



(I) 



(3) 



where B' is the magnitude of the field outside the sheet, 5' is 
a measure of the sheet thickness and the field is in the direc- 
tion parallel to the sheet and perpendicular to the direction of 
current flow, which is locally considered to be Z'. The particle 
population in the sheet consists of two counter-drifting relativis- 
tic Maxwellian distributions the density of which falls with X' , 
and whose drift provides the net current. From the pressure bal- 
ance across the sheet f ollows the dependence of the density of 
each distribution on X' (Kirk & Skiaeraasen 2003): 



NL = 



N' ±0 cosh 



±0 



(4) 



(5) 



\6nmc 2 ® 

where = k^T /(mc 2 ) is a dimensionless temperature asso- 
ciated with the particle distribution. The condition connecting 
the particle density, temperature and current sheet thickness is 
dKirk & Skiaraasedl2003HLvubarskv& Kirk 200jl): 



@mc 



4nN' ±0 e 2 y ± 



= 6"Bl 



(6) 



where y ± and B ± are the Lorentz factor and the corresponding 
speed, normalized to the speed of light, of the drift of the dis- 
tributions, in the wind frame. Finally, by assuming that the ideal 
gas law holds for the distribution we get: 



2N' ± T' 



(f - l)(e' - 2N' ± mc 2 ) 



(7) 



where e' is the energy density associated with the two counter- 
streaming distributions. For the relativistic Harris sheet solution, 
this equation becomes: 



r 



l + 







l (i +£) +r ± (Ki(im/K 2 (im - 1) 



(8) 



Y 



where K\ and K2 are modified Bessel functions of the second 
kind. In the inner part of the pulsar wind, which is of interest 
to us, the thermal particles in the distribution are predicted to 
be highly relativistic, something which allows us to set y = 4/3 
in the above equation and also has as a result that the thermal 
energy of the particles in the distributions will greatly exceed 
their rest-mass energy: » 1. Using the approximations 
-1 



..-2 



Ki{x<& 1) 
K 2 (x «1) 2x~ 

Eq. |8]can be simplified to: 
1 

■ 



(9) 
(10) 
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This means that the drift of the distributions in the current sheet 
is not relativistic, and beaming effects caused by it can be ig- 
nored when calculating the emitted radiation. Therefore the dis- 
tributions can be for this purpose considered to be isotropic in 
the frame in which the sheet is at rest. Using this approximation 
we can solve Eqs.[5]and[6]for © and N' ±Q . In doing this, we have 
to take into account the fact that the angle arccos(n ■ r) between 
the sheet normal and the radius unit vector changes when mov- 
ing from the lab to the wind frame. 

The magnetic field from the solution of lBogovalovl (1 19991) is 
always parallel to the sheet, if one ignores reconnection effects. 
This means that close to the light cylinder and for a radial out- 
flow one cannot ignore the radial component of the field. The 
full field in the wind frame will be: 



(12) 



and the toroidal component will prevail in this frame only 
for R > F (we will call this region the "far wind" region). 
Conversely, in the region R < F (which we will call the "near 
wind" region), it is a good approximation to ignore the second 
term under the square root and approximate B' Q B^c/R 2 - 

If one considers a perpendicular rotator (a pulsar with x = 
7r/2), the temperature and density in a current sheet segment are 
given by the simple expressions: 




= 



/ «lcA sin ft 
\ 2R 

\6nmc 1 ® 



2/3 



(13) 
(14) 



The dimensionless parameter «lc appearing above is called the 
strength parameter, and is defined at the pulsar's light cylinder 

as: 



eB hC P 
2nmc 



(15) 



Taking a pulsar's moment of inertia to be / = 10 gem , the 
spindown lu minosity is calcula ted by the period and the period 
derivative as (lAbdo et al.ll2010t) : 



£33 = 4tt z 10 ii PP- 



(16) 



where we have normalized the luminosity to the value 10 33 erg/s. 
The magnetic field (in Gauss units) can then be expressed as: 



B LC = 46.83£ 3 1 f P 



The strength parameter then can be expressed as a function of 
the spindown luminosity only: 



a LC = 1.3 x lO 6 ^ 



8*1/2 



33 



(18) 



All GeV pulsars detected so far have strength parameters in the 
range 10 s - 10 u . 

Finally, A < 1 is the fraction of a half wavelength that a sheet 
occupies, measured in the radial direction (i.e. not perpendicu- 
larly to the local sheet plane), in the lab frame. We will assume 
that A is constant with radius for the sake of simplicity, how- 
ever it is far from clear how this parameter evolves with radius 
and obliquity y close to the light cylinder. We will also assume 
that F is constant, although generally reconnection in the cur- 
rent sheet has been shown to increase T with radius in the fa r 
wind zone dLvubarskv & KMdl2001t Ikirk & Skiaeraasenll2003l) . 



a result that might also apply to the near wind zone. These as- 
sumptions are based on the fact that only a very limited radius 
interval contributes to the gamma-ray radiation, as we will see 
below, therefore if the change in A and T happens with a scale 
larger than r-uc, it is not relevant to the present estimations. 
Substituting for «lc in Eq.[13]we get: 







1.6 x 10 5 £^ 3 



A sin ft 
R 



2/3 



(19) 



The requirement that » 1 in the near wind region R < Y 
translates to: 



A» 10 



-8 



E\j 2 sin ft 



(20) 



which, as we will see in the following, is easily satisfied for all 
y-ray pulsars, provided the emitting region is not very close to 
the polar axis. In any case, the validity of the » 1 assumption 
has to be checked a posteriori for all cases of studied pulsars. 

In the following, all results will refer to the near wind region, 
R < T, but still outside the pulsar's light cylinder, close to which 
we assume that the current sheet is formed. Therefore we will 
consider a minimum value for the radius of the radiating sheet 
of Rmm = 1 ■ Our model is valid only for relativistic outflows, so 
we will also set a lower limit on T of r m ; n = 10. 

In the next paragraphs, we assume a perpendicular rotator 
(X — 7i IT) in order to give some analytical estimates concerning 
the emitted radiation and the characteristics of the current sheet 
and its populations. These estimates can be generalized to the 
oblique rotator (see Appendix |B), however the conclusions re- 
main essentially the same as long as the line of sight is not close 
to the edge of the current sheet, where £ = tt/2-X- 

2.2. Larmor radius of hot particles 

An important condition for the consistency of our model comes 
from the requirement that the hot particles' Larmor radius in the 
full field between the current sheets should be smaller than the 
sheet width. In the wind frame this requirement translates to: 



(y)mc 2 

~7bY 



» Ar LC R 



-1/2 



(21) 



This inequality is equivalent to 0'^ 2 » 1 which should hold 
close to the light cylinder, by the assumptions of the model. 



(17) 2.3. Energy of emitted photons 



The mean Lorentz factor of the electron/positron distribution in 
the current sheet is (y) ~ 30. Electrons of this (y) gyrating in 
the field B' () , radiate photons of energy (in the observer's frame): 



3 , eB' 

s * -r(30) 2 fc— 

2 mc 



(22) 



After some manipulation, one finds the highest observable pho- 
ton energy: 



->GeV 



2X 10 



_ 4 El' 6 r(Asin0 4/3 



R 



10/3 



(23) 



where the photon energy is measured in GeV. In Eq.|23]we have 
kept only the poloidal component of the magnetic field, since it 
dominates close to the light cylinder, at R < T. All pulsars in the 
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first Fermi/LAT catalogue have Eli IP » 1- This is why it is 
possible for most objects in the catalogue to find an appropriate 
combination of F, A and R m \ n > 1 which will bring ^Gev to the 
regime observed by LAT, while at the same time satisfying all the 
restrictions on A, T and R. This is especially easy for the highest 
luminosity pulsars (£33 » 1), or for the very low period ones 
(P « 1 s, millisecond pulsars), as one can deduce by inspecting 
Eq. [23] In other words, for pulsars of the same spindown lumi- 
nosity, millisecond pulsars are more likely to reach GeV ener- 
gies. Alternatively, if T and A are the same, millisecond pulsars 
with lower spindown luminosities than young pulsars can emit 
in the FermifLPS band, a trend that might be observable as our 
statistics on gamma-ray pulsars increase. 

From the constraints on A and R, we can calculate an abso- 
lute maximum on the peak energy of the emitted spectrum, for 
$ = nil: 



^GeV.max - 2.3 X 10 T 



^33 
P 



it. The pressure in the sheet will fall, thus causing its com- 
pression by the external magnetic field. This can initiate com- 
pression driven magnetic reconnection, a phenomenon pre- 
viously studied in the context of the interaction of current 
sheets in a pulsar wind with the wind's termination shock 



(lLvubarskvl2003l:lPetri & Lvubarsky|2007l:lLvubarskv & Livertsl 



(24) 



120081: ISironi & Spitkovskvl201 1). Reconnection has been shown 
to be able to accelerate particles to non- thermal distributions 
above the thermal peak of the particle spectrum (and also, as al- 
ready mentioned, to accelerate the bulk flow and cause the sheet 
width to rise). A detailed study of this radiation loss induced pro- 
cess is beyond the scope of this work. It should be stressed, how- 
ever, that if t[/t' R < 1, some acceleration processes ought to be at 
work in order to supply the sheet with energetic particles. This 
will in all likelihood result in power-law distributions in the cur- 
rent sheet, a signature that might be observable in the spectrum 
above the MeV-GeV peak for the highest luminosity pulsars and 
for those millisecond pulsars with t' s /t' R < 1 (something that has 
not been modelled in the present article). 



As we will see below, A is often more severely constrained by 
the radiation reaction limit, so the energy £GeV,max will not nec- 
essarily be reached for all objects. 

The rapid fall of the peak photon energy with radius means 
that in this simple model the most energetic photons come from 
close to the light cylinder. In a realistic situation it is to be ex- 
pected that either F or A or both will rise with radius, making 
the fall of fiGeV with radius less steep, however even if A and F 
proved to rise linearly with R, the emitted photon energy would 
still fall as oc R l . 



2.4. Energy losses 

For the relativistic Harris equilibrium to hold in a quasi-steady 
state, equilibrium has to be established in the wind frame within 
a timescale of the order of magnitude t' R ~ R(ra>y l , which is the 
timescale on which the magnetic field change is comparable to 
its magnitude 6B' ~ B' . Equilibrium is established in the current 
sheet within a timescale 



Ar LC 



4,5 



Therefore the requirement f' < t' R holds close to the light cylin- 
der as long as A < r 1 , while this condition is relaxed with radius 
to A < R/F. 

For the particles not to suffer catastrophic energy losses, 
their synchrotron cooling timescale t' s should be larger than t' R . 
Comparing the two one finds: 



4 4.57 



TR n/3 



£4/3 (Asin<D 2/3 



(25) 



For the lower £33 pulsars the inequality t' s > generally holds 
for R > R m i n , but this may not be the case for some millisecond 
pulsars, because of the linear dependency of the ratio on pul- 
sar period P. The steep radius dependence, however, relaxes this 
condition rapidly, so that when it breaks down, it does so only 
for a very limited R -range close to the light cylinder. 

If t' s /t' R < 1, as is often the case for high spindown pul- 
sars and millisecond pulsars, particles will lose energy rapidly 
and cool, disturbing the pressure equilibrium between the pop- 
ulations in the current sheet and the magnetic field outside 



2.5. Radiation reaction limit and maximum emitted energy 

If the particles of mean energy ~ 3&mc 2 in the current are gain- 
ing energy through acceleration in an electric field E ~ %B with 
£ < 1, then their acceleration is limited by radiation losses. For 
radiation-reaction limited synchrotron emission (in the full field 
amplitude B' ), energy losses compensate possible energy gain: 

dW 

— = e£B' c < eB' c (26) 

Jsyn 

From this expression one can calculate an upper limit for A: 

fl5/2p3/4 



A < Aiim = 256 



El{ & sin ( 



(27) 



'33 



The £33 dependence in the denominator implies that for the most 
powerful pulsars, the sheet is thinner close to the light cylinder. 
For the Crab we get Ai; m = 3 x 10~ 4 sin £ at R m i n , whereas 
for most of the weaker pulsars with £33 < 10 2 , Aij m is greater 
than unity, therefore not presenting a constraint. All millisecond 
pulsars detected by £erm//LAT fall into this category. 

Inserting Aij m in Eq. [23] one recovers the classical syn- 
chrotron limit: 



eV .max 



< 0.36F 



(28) 



where the factor F comes from the boosting of the photon's en- 
ergy to the lab frame. The minimum of Eqs. 1241281 then, can 
give an estimate of the maximum value of the GeV cutoff for 
any individual object. 

Inserting Au m into Eq.[25] one gets a lower limit for the ratio 
t' s /t R (applicable mainly to higher-spindown pulsars, for which 
Aiim < 1): 



0.23rR 2 P l/2 E^' 4 



(29) 



We see that, mainly for young pulsars of large spindown lumi- 
nosity, the right hand side can be less than unity. These pulsars 
either have smaller A in order to keep t' s /t' R > 1, or rapid reac- 
celeration mechanisms in their current sheets as argued above 
(or both). Therefore there should be a trend in the P — P dia- 
gram to observe more objects with acceleration signatures (i.e. 
power-law tails) in their spectrum as one moves to higher spin- 
down luminosities, or equivalently as one moves to higher P and 
lower P (upward left region in a P - P diagram). 
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Fig. 1. Sky maps of the pulsar's total luminosity shown for two different bulk Lorentz factors T = 10 and F = 50 and for different 
X- In the upper panels x — 7r /2, in the middle ones x — 7r /3 and in the lower ones x — tt/6. All the plots have been made using the 
parameters A = 0.01 and R m \ n — 1. The luminosity is given by integrating £ v (Eq. IA.4b over all frequencies and multiplying with 
d 1 , where d is the distance of the pulsar from the observer. A beaming factor of To. = nas been used. 



3. Lightcurves, spectra and y-ray luminosity 

3.1. Lightcurves 

The expected lightcurves and spectra produced in the current 
sheets of a pulsar can be computed numerically by integrating 
the emission coefficient of the particles in the magnetic field of 
the current sheet along the line of sight to the observer. The pro- 
cedure is explained in the Appendix. 



In the wind model of the pulsar radiation two pulses per pe- 
riod are expected, the separation and the width of which varies 
with the obliquity x and the line of sight angle (. The width of 
the pulses depends on the bulk speed of the outfl ow, with wider 
pulses for lower F (Kir k et alj|2002t lPetrill20ill) . Examples of 
the variations of lightcurves with^, £ and T are shown in Fig.Q] 
where we have used a model pulsar of surface magnetic field 
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equal to B — 10 G and period P = 0.01 s in order to plot the 
sky maps of the pulsed total luminosity. 

In the left column of Fig. [TJ sky maps are shown for F = 10, 
while in the right column the value of the bulk Lorentz factor 
is F = 50. Fhe width of the pulse depends strongly on F, as 
expected. 

Fhe upper maps in both columns represent the case of the 
perpendicular rotator, where x — 7r /2, the middle ones are for 
X = tt/3 and the lower ones for^ = n/6. As the obliquity falls, 
the width of the pulse increases while the peak luminosity falls, 
resulting in the overall phase-averaged luminosity remaining es- 
sentially constant for different obliquities. Fhe widest pulses are 
observed for the smallest inclination angles. However, for very 
small x the line of sight to the observer has to lie very close to 
the equatorial plane of the pulsar for the pulse to be observed, 
since the luminosity is significant only for f < n/2 — x- This is a 
likely geometry for the objects which are observed to have two 
wide pulses with a phase separation of dcf> ~ 0.5. 

From the sky maps in the case of F = 10 it can be discerned 
that there is a slight substructure in the pulses. There is a slight 
dip in luminosity at the peak of each pulse, resulting in two sub- 
peaks appearing for all <f. This is caused by the structure of the 
current sheet. In the middle of the sheet the magnetic field is 
exactly zero, rising towards the sides, while the density of the 
hot particles is highest and falls towards the sides. Therefore the 
bulk of the radiation of each current sheet segment comes from 
two regions away from the sheet midplane, where the product of 
particle density and magnetic field is largest. The double-peaked 
shape of the pulse reflects the two luminosity peaks within the 
sheet. Because of the shrinking pulse thickness with rising F 
this phenomenon should be observable only for the lowest bulk 
Lorentz factors. 

It is useful to note here that the highest luminosities in the 
sky maps are correlated with the highest peak energies £GeV, 
resulting in pulses that get sharper for higher energy photons. 
This is shown in Fig. [2j where the flux vF v has been plotted in a 
color map as a function of emitted energy & (in eV) and phase <p. 
A horizontal slice of the map shows the lightcurve at a specific 
energy while a vertical slice shows the spectrum at a given phase. 
The pulses become narrower with increasing energy, as can be 
seen from the wider flux variation in a horizontal slice as one 
moves higher in the map. Also, the peak energy varies by many 
orders of magnitude within one phase, and is lowest between 
pulses, something that can be seen when comparing spectra at 
different phases. 

In Fig. [3] we show how the lightcurves change when moving 
the minimum radius R m i n from the light cylinder to a distance in 
the "far wind" region, R m \ n > F. For low R mm , close to the light 
cylinder, the pulse shape is almost symmetrical with respect to 
its peak (also seen in Fig. [TJ. As one moves towards larger R m i n , 
the pulses stop being symmetrical and instead present an abrupt 
rise in luminosity followed by a gradual fall, making them asym- 
metrical with respect to the peak. Also, the total luminosity in 
the pulse is reduced by many orders of magnitude as one moves 
outwards in the wind. In the example of Fig. [3] a change from 
^min = 1 to R m \ n = 30 results in the decrease in luminosity by 
seven orders of magnitude. These results agree with previous 
investigations of pulses from the far wind region (for example 
iKirket all (120021) ). 

From the sky maps we have shown we can conclude that the 
separation of the two pulses depends on x an d f> with a sep- 
aration of A(p ~ 0.5 when the line of sight lies at the equato- 
rial plane. However, as we have pointed out above, the phase- 
averaged luminosity changes significantly only if the line of 



vF v d 2 (erg/sec) 




0.2 0.4 0.6 0.8 1 

♦ 



Fig. 2. The peak flux v/ r v (multiplied by the square of the distance 
to give a luminosity estimate) as a function of phase <p and peak 
energy S. S ranges from 10 MeV to 1 TeV. The same pulsar 
parameters were used as in Fig.[TJ with^- = 7r/3 and £ = tt/2. 



sight lies above the wind, i.e. if f < 7r/2 — %. In the opposite 
case the overall luminosity is not sensitive to variations of^- and 

3.2. Phase-averaged spectra 

In Fig.[4]we have used the same model pulsar as in the previous 
section to calculate phase averaged spectra and their dependence 
on the model parameters. In the upper plot we show the varia- 
tion of the spectrum when changing the bulk Lorentz factor F. 
As is expected from Eq. [23] the peak of the spectrum rises lin- 
early with T, however the overall luminosity falls with rising T. 
This can be attributed to the diminishing of the area of the sheet 
from which doppler-boosted radiation is received: the boosted 
area scales as 1 /F 2 , something that is only partly compensated 
for by the square of the doppler factor in the calculation of the re- 
ceived flux (given in the Appendix). The most strongly boosted 
radiation comes from the line of sight, however the luminosity of 
the region directly intersecting the line of sight is not the major 
contribution to the overall observed luminosity. This is caused 
by the fact that one observes only the effects of the perpendicu- 
lar (i.e. azimuthal) field at # = £, while the effects of the much 
larger poloidal field B r come from a region at angle ~ 1 /T to the 
line of sight. Therefore when the line of sight is directly aligned 
to the edge of the sheet, where £ = n/2 - x, still two pulses 
of the same amplitude are observed, coming from the two parts 
of the folded current sheet that are at angle ~ 1 /F to the line 
of sight, whereas the edge itself contributes only a little to the 
over all flux. Thi s is different to what has been predicted in the 
past (lPetrill2011h . However, a single wide pulse will be observed 
when looking over the edge of the sheet, the disadvantage in this 
case being that the luminosity decreases quickly with decreasing 
f making such pulses difficult to detect. 

In the middle panel the change in the phase-averaged spec- 
trum is shown when varying A. The dramatic rise of the cutoff 
energy and of the luminosity with A can be explained by the de- 
pendency of the peak of the spectrum on A, given in Eq. [23] as 
well as the rise of the relativistic temperature in the current sheet 
with A as seen in Eq. [19] 

In the lower panel we demonstrate the dependence of the 
cutoff and luminosity on the minimum radius R m \ n . The peak fre- 
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JC = n/2, ?=ti/2, A= 0.01, R min =1 
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Fig. 3. Variation of the pulsar's lightcurves for different R m i n . 
The bulk Lorentz factor in these examples is F = 10 and the 
obliquity is^ = n/3. The rest of the parameters are the same as 
in Fig.Q] 



quency of the spectrum and the overall luminosity both rise very 
strongly with decreasing R m i n , something that is in agreement 
with Eq. [23] and is also expected because of the dependence of 
the temperature of the sheet particles on radius, given in Eq. [T9j 
Therefore the bulk of the received high-energy radiation comes 
from a region of limited range in R close to ^ m j n . 
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Fig. 4. Phase-averaged spectra for the same model pulsar that 
was used in Figs. Q] and [3] All the plots have been made using 
X = n/2 (perpendicular rotator) and £ = 7r/2. Here the variation 
of the phase-average flux with the parameters F, A and R^ is 
shown. 



3.3. Examples 

In Fig. [5] we present two examples of spectra and lightcurves 
of Fermi/LAiT detected pulsars, calculated using our model. 
The first example is of the millisecond pulsar PSR J1614-2230, 
which has a period of P — 3.2ms and a spindown luminosity of 
£33 = 5. Its spectrum calculated according to our model can be 
seen in the upper left panel of Fig. [5] alo ng with the best fit p ower 
law with exponential cutoff, as given in lAbdo et al. I J2009I) . The 
parameters used were A = 0.25, F = 20, x - 7r /2 and ( = n/2. 
The lightcurve corresponding to the same parameters is seen in 
the lower left panel. The phase of the two peaks is correctly pre- 
dicted, however the width of the simulated pulses is narrower 
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than what is observed. This is because in this example A > 1 /F, 
and in this case the lightcurve cannot be accurately predicted by 
our model (as explained in Appendix iBt. 

The second example, shown in the plots on the right of 
Fig. [5] is of pulsar PSR J0205+6449, which was discovered 
by Fermi/LAT in the supernova remnant CTA1. This pulsar has 
a high spindown luminosity of £33 = 2.7 x 10 4 and a pe- 
riod of P — 65.7ms. It is a young, energetic pulsar for which 
Aiim ^ 1.5 x 1CT 3 . The parameters that we have used for these 
plots are A = Aii m , T = 80,^ = 7r/5 and £ = tt/2. The relatively 
high Lorentz factor needed to reach the cutoff of £ » 3GeV 
causes the pulses to be rather sharp, something that is also ob- 
served with Fermi/LAT. 

From these two examples as well as from the phase-averaged 
spectra shown in Fig . |4] it can be seen that our model predicts a 
less steep fall and a substantially larger energy flux at photon 
energies above 10 GeV than would be expected using a power- 
law plus exponential cutoff fit, something that implies that for 
the more energetic pulsars, which have larger bulk Lorentz fac- 
tors, one could detect radiation in the sub-TeV to TeV regime 
explored by Cerenkov telescope arrays. Specifically a good fit 
to our spectra is given by a power-law with a sub-exponential 
cutoff: 



dN r-P 
— oc & 1 exp 

do 



-•cutoff 



(30) 



Our model predicts a cutoff with b = 0.35. This index is not 
sensitive to the obliquity x or the observer angle as long as 
the line of sight intersects the equatorial current sheet, i.e. when 
£ > 7t/2~x- The index b can reach the value ~ 0.4 for ( < n/2~x, 
however the very low luminosities that are expected in that case 
make the majority of those objects unobservable. 



3.4. Peak energies and y-ray luminosities 

In Fig. [6] we show a prediction of the region in the P — P where 
spectra that peak in the range 0. 1 - 10 GeV are espected, accord- 
ing to our model. Two shaded regions are shown that correspond 
to^A = 0.1 and T = 100 (yellow) and A = 0.001 and T = 100 
(green). We have also plotted the lower and upper borders of the 
region that corresponds to A = 0.1 and T = 10 (red lines). It is 
obvious that almost all pulsars already detected by Fermi/LAT 
fall in at least one of these regions, leading to the conclusion 
that one can accommodate the cutoffs of almost all pulsars using 
parameters in the range 0.001 < A < 0.1 and 10 < F < 100. 

In principle the peak energies emitted by the near wind re- 
gion of a pulsar wind are not constrained to the GeV regime. 
Peaks at TeV energies are likely excluded by the synchrotron 
limit given in Eq. [28] however there should be a population 
of pulsars with spectra peaking in the MeV regime (or lower). 
These pulsars can be found in a region of the P - P diagram 
roughly below the yellow shaded region of the P — P diagram in 
Fig. [6] The population of MeV pulsars will include objects with 
spindown luminosities lower than the ones which are detected 
by Fermi/LAT. This becomes obvious by inspecting Eq. [23] for 
lower spindown luminosities the same T and A lead to lower 
peak energies. For example a pulsar which is a perpendicular 
rotator with P = 0.5 s and P = 10~ 15 with bulk Lorentz factor 
F = 10 and A = 0. 1 will radiate a spectrum that peaks at energies 
8 =* 0.3 MeV. 



The overall luminosity emitted by the near region of a pulsar 
wind from a perpendicular rotator according to our model can be 
estimated to within a factor of two by the expression: 



-£tot — 



2.7 x 10 



31 



(t^t 1 (JLt 71 (A. 
I10 3 / lo.i; 0.1 



(31) 



where Z?lc is given in G and P in s. In terms of spindown lumi- 
nosity the above can be expressed as: 



7.6 x 10 



-0.9 



£33' I „ . 



A 

oTi 



(32) 



It is obvious from the above expression that, since no trend is 
observed for y-ray pulsars to have a y-ray luminosity that is anti- 
correlated to their spindown luminosity, the values of T and A 
should vary significantly from object to object. 

If the peak of the spectrum falls into the y-ray regime, then 
Eq.(3jJ gives a good estimate of the y-ray luminosity of the pul- 
sar. The dependence of _£tot on% and £ is very weak as long as the 
line of sight intersects the wind, £ > n/2 — x, therefore one can 
use the above formula to reach approximate conclusions about 
the y-ray luminosity of a pulsar. Eq s . [3T1 and 1231 (which gives us 
an estimate of the peak energy of the phase-averaged spectrum) 
show that, once P and P are known, the emitted luminosity and 
peak energy can be fixed essentially by the two parameters A and 
r. These two equations can also be inverted to deduce values of 
A and F for a sample of pulsars for which P, P, X y and fiGeV 
are known. This was done for the pulsars of the First Fermi/LAT 
pulsar catalog and the results can be seen in Fig. [7] where the A 
and T have been plotted for each object of the sample. 

The blue squares represent the young pulsars in the catalog 
while the red triangles are the millisecond pulsars. For most ob- 
jects A ranges between 10~ 3 and 1, and a trend is observed for the 
millisecond pulsars to have larger values of A, slightly smaller 
than unity. The bulk Lorentz factor F ranges roughly between 10 
and 100, with most objects having relatively low Lorentz fac- 
tors in the range 10-30. For the most part the deduced T and A 
are reasonable, except for a region where values A > 1 are pre- 
dicted, which are unphysical since A has to be less than unity 
for the current sheet to be narrower than the wind's half wave- 
length. This is a result of the fact that the y-ray luminosity of 
these objects is relatively high, so that by Eq. (3jJ a higher A is 
needed. Our model is possibly too simple to accomodate all of 
the pulsars in the catalog, and maybe there is a contribution from 
another physical mechanism or emission region which results in 
higher luminosities in the exceptional cases of the objects with 
A > 1. 

In the diagram on the right of Fig. [6]we have plotted the y- 
ray luminosity L y of two random samples of pulsars (millisecond 
and young pulsars) with respect to their spindown luminosity 
L sp . To this purpose two random samples of pulsars were used, 
with periods and period derivatives in the ranges -2 < log P < 
and -15 < logP < -11 (for the young pulsars) and -3 < 
logP < -2and -20 < logP < -17 (forthe millisecond pulsars). 
We have chosen £ and x randomly in the interval [0,7r/2], the 
bulk Lorentz factor in the interval 1 < log T < 2 and the sheet 
thickness either within log An m - 2 < log A < log Ai; m in the case 
Aiim < 0.1 or within -3 < log A < -1 if Ai im > 0.1. The y-ray 
luminosity was calculated by integrating the phase-averaged flux 
between 100 MeV and 100 GeV and using a beam correction 
factor 'F n - 1, in order to have a direct comparison with the 



results of lAbdo et alJ (120101) . 

About 12% of all young pulsars in our sample and 30% of 
the millisecond pulsars have a cutoff within the inteval 100 MeV 
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Fig. 5. The phase-averaged spectrum and simulated light rcurve of millisec ond pulsar J1614-2230 as calculated using our model, 
along with the best power-law plus exponential cutoff fit dAbdo et al.ll2009l) . The parameters that were used are A = 0.25, F = 20, 
X = £ = x/2 for the millisecond pulsar (left). For the young pulsar in the CTA 1 supernova remnant we used A = Au m , T = 80, 
X = n/5 and £ = n/2. 



- 10 GeV and are shown in Fig. [6] These percentages reflect our 
choice of the range in A. We have constrained A to be within two 
order of magnitudes of the value min(A m j n , 0.1). If the range of A 
were narrower, resulting in systematically larger values, then the 
percentages quoted above would rise. Lower choices of A lead 
to lower £GeV so there are more pulsars that have peaks below 
100 MeV and therefore don't make the cut. 



The large beam correction factor used for this plot has as a 
result the occurrence in some rare cases of objects with L y > L sp , 
whereas in most cases the y-ray luminosity is a few percent of 
L sp . The young pulsars tend to lie higher in the diagram than mil- 
lisecond pulsars, something that reflects their larger spindown 
luminosity in combination with the fact that the gamma-ray effi- 
ciency is similar in both samples. 

In lAbdo etalJ i ]2010l) a trend was reported that the gamma- 
ray luminosity is proportional to the square root of the spin- 

1 II 

down luminosity for the highest luminosities: L y oc L s i . We 
don't have a clear indication of such an effect in our simulated 
sample, meaning that the gamma-ray efficiency of the highest 
spindown pulsars is predicted to be similar to the one of the 
lower spindown ones. The best fit line to the points in Fig. [6] 
gives L y oc Ls p 95 for both young and millisecond pulsars and if 
the statistical error is taken into account, then the relationship is 
compatible with the proportionality L y oc L sp . 



4. Discussion and conclusions 

In this paper we have shown how gamma-ray pulses can natu- 
rally arise within the framework of the pulsar wind's equatorial 
current sheet outside (but close to) the light cylinder. The advan- 
tage of this emission model is that it is an intrinsic mechanism 
that naturally produces peak energies in the MeV-GeV range. It 
can give meaningful results for almost all pulsars, irrelevant of 
age or environment, employing only a few parameters: the bulk 
Lorentz factor of the outflow T, the sheet thickness A, the obliq- 
uity x an d the angle of the rotational axis to the line of sight £. 
To these the minimum radius R m i n can be added (which, in the 
calculated examples, has been set equal to ^mi n = 1). However 
it is important to note that R m { n has to be sufficiently close to the 
light cylinder, so that the peak energy of the spectrum reaches 
the GeV regime. The parameters of the model are constrained 
by the characteristics of the pulsar, and by the physics of the 
current sheet. The novelty of our model is that it makes an at- 
tempt to account for the previously ignored region of the wind 
between R = 1 and R T, and is therefore able to take advantage 
of the strong poloidal field close to the light cylinder, something 
that has not been discussed in previous wind models. 
The predictions of our model are the following: 

1. Two pulses per pulsar period are expected. The pulses have 
the same amplitude. 

2. The shape of the pulses is symmetric with respect to their 
peak, as long as emission starts close to the light cylinder, 
and become increasingly asymmetric with rising R m i n - 



9 



I. Arka & G. Dubus: Pulsed y-rays from thermal particles in pulsar wind current sheets 



0.1<E GeV <10 




39 
38 
37 
36 
35 
34 
33 
32 
31 
30 



33 



Gamma-ray luminosity vs spindown luminosity 





















■ 


>••""""' 

A 












& 

»*' A- 




^A ■ 






A itll ,. 

„ "" / 


A A 
\ A ■ % " - ■ 

f^^At A 


■/■^"m- a 


■ 




(i """"a 

A 4 
A a aA 


" A V 

A k A .4 


AA^ 

A 




'"I*..-- - " 




*h -* A 1 


A ■ 
■ A 










: a'*t a - 

A A 








A 












A 













34 



35 



36 

Log 10 (L sp )(erg/sec) 



37 



38 



39 



young pulsars 
L v =L cri 



Fig. 6. On the left: a P - P diagram including the pulsars observed by Fermi/LAT and a plot of the predicted y-ray luminosity in the 
range 0.1-10 GeV for a random sample of pulsars. On the right: the gamma-ray luminosity L y versus the spindown luminosity L sp 
for two random samples of millisecond and young pulsars. The lines L y 
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Fig. 7. The parameters A and F estimated for the pulsars in the 
FermifLAT one-year catalog. Blue squares are the young pulsars 
while red triangles represent the millisecond pulsars. 

3. There is no significant interpulse emission in the GeV range. 

4. The width of the pulse decreases with increasing photon en- 
ergy. 

5. The separation of the peaks varies with the obliquity x an d 
the angle to the line of sight as in previous models of 
pulsed radiation from the pulsar wind. 

The features of the pulsed emission that cannot be explained 
by our model, such as the interpulse or the peaks of different 



intensity, could be accomodated by a model in which radiation 
comes both from within the light cylinder, possibly in an outer 
gap, as well as from the beginning of the wind. The difficulty of 
ou ter gap models to produce double peaked lightcurves, as noted 
in iBai & Spitkovskvl (l2010h . in combination with the difficulty 
of our near wind model to produce single pulses, point to the 
need for a combined model in which the high energy emission 
starts within the light cylinder and continues in the equatorial 
current sheet outside the light cylinder. Therefore the physics of 
the region of formation of the current sheet is very important in 
the endeavour to understand y-ray emission from pulsars. 

For pulsars in binaries, the possibility exists that the favoured 
radiation process is not synchrotron but rather inverse Compton 
o n the low frequency ph otons of the massive star, as investigated 
in iPetri &Dubusl(l201ll) . In this case the peak energies and emit- 
ted luminosities by the two mechanisms can be comparable, de- 
pending on whether the comoving energy density of the low fre- 
quency photons in the wind frame is comparable to the magnetic 
field energy density close to the light cylinder. For the pulsar 
B 1259-63, which is a member of a binary with a B2e star, syn- 
chrotron radiation close to the light cylinder would give a peak 
energy £GeV ~ 3 (assuming T — 10 and A = 0.01), while inverse 
Compton radiation on the star's photon gives a peak energy that 
is constrained by the Klein-Nishina limit to £GeV ~ 10. In this 
case the low frequency photon energy density is several orders 
of magnitude lower than the one of the magnetic field, therefore 
synchrotron radiation should be the dominant component. This, 
however, is not necessarily the case for other binary systems, 
such as LS 5039, where the seed photon energy density exceeds 
the one in PSR B 1259-63 by a factor of 200. The detection of 
orbital modulation of the received y-ray flux can be the criterion 
on which to distinguish between the two emission mechanisms. 

A po ssible issue with o ur model is that the supersonic so- 
lution of iBogovalovl (1 19991) might not apply for the objects of 
lower bulk Lorentz factors. Particularly, the wind can acceler- 
ate, as is pred icted in various magnetohydrodynam ic models of 
pulsar winds dBeskin et al.lll998t Kirk et alJl2009h . In this case 
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the shape of the current sheet can in reality be different, since 
there is generally also a polar magnetic field component, which 
in our model was zero. However, as was mentioned above, the 
very steep dependence of £GeV on R implies that even if the 
parameters A and F changed with R, the present model would 
still be able to explain the main features of the observed spectra. 
Adding to these considerations the fact that the particle density 
within the current sheet and the magnetic field decrease as the 
wind expands to larger R, the conclusion is reached that the emit- 
ted luminosity, which depends on N' and B' Q will also decrease 
with radius, and this in combination with the decreasing peak of 
the spectrum means that the main contribution to the spectrum 
near the cutoff will still come from a very limited region close 
to R m \ n . Therefore, possible changes in the dynamic evolution of 
the wind should not have a large effect to the observed spectrum 
and luminosity, while the exact shape of the current sheet should 
influence mainly the shape of the pulse, but not the cutoff energy 
or the emitted luminosity. 

It is expected that reconnection will generally alter the 
physics of the current sheet beyond the pulsar's light cylinder, 
and might result in an evolution of the current sheet thickness 
and the wind Lorentz factor, especially in objects for which 
the radiative timescales of the thermal particles are very short 
in comparison to the evolution timescale of the magnetic field. 
Furthermore, particle acceleration takes place during the recon- 
nection process which might result in a non-thermal tail to the 
thermal distribution, the characteristics of which will change as 
the current sheet evolves with radius. These phenomena will 
leave their imprint on the pulsed high-energy spectrum, and 
need to be investigated in a self-consistent way, something that 
is beyond the scope of this article. We can, however, give a 
simple order-of-magnitude argument about the reconnection- 
accelerated particles: deep in the current sheet particles can be 
accelerated by reconnection electric fields to energies higher 
than the thermal peak. These particles would have Lorentz fac- 
tors extending to y ~ ^a^c/F (in the wind frame), thus giving 
rise to photons up to an energy 

£cev~8.8-^ (33) 

where f = E'/B' < 1. This could extend beyond the thermal 
peak for large enough £ or £33, whereas millisecond pulsars are 
again favoured by the dependence on the inverse of the period. If 
these high energy non-thermal particles escape the acceleration 
site and radiate in the field within the current sheet, their emis- 
sion could give rise to a power-law t ail extending beyond the 
GeV cutoff dZenitani & Hoshinoll2008l) . This mechanism should 
be prominent for objects for which t' s /t' R < 1, which tend to clus- 
ter at the upper left part of the P — P diagram, as discussed. It is 
for these pulsars that acceleration by reconnection could become 
prominent and be observed in the form of power-law tails in the 
GeV-TeV regime. 

Another possibility is that the particles in the current sheet 
are already accelerated to a power-law distribution when the 
sheet starts radiating, in which case a different distribution 
function would h ave to be used in order to c alculate the current 
sheet parameters dBaUkhin & Gedalir]l2008l) . This might apply 
particularly to millisecond pulsars, which have lower surface 
magnetic fields which lead to lower pair production rates, and 
therefore less dense plasmas in their magnetospheres. In this 
case a non-thermal particle distribution seems likely to describe 
the physics of the current sheet in a more consistent way. We 
defer the investigation of such particle distributions to a future 



article. 
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Appendix A: Calculation of received synchrotron 
spectrum by a relativistcally moving sheet 

The flux that an observer receives fr om a relativistically mov ing 
source is calculated by the formula dLind & Blandfordlli985h : 




(A.l) 



- Primed quantities are in the wind frame, i.e. the frame where 
the outflow is at rest. Since the wind is assumed to be flow- 
ing radially, at each azimuth and polar angle (<p, #) there is 
a different, local, radially moving "wind frame". Quantities 
pertaining to a point of the wind with coordinates (r, ft, <p) are 
calculated in this local wind frame. 

- When calculating the emission coefficient / we should take 
care to calculate it in the direction of the line of sight, taking 
into account the aberration of photons in the wind frame (we 
will return to this matter later). 

- The volume element is taken in the observer's frame (or 
the frame in which the pulsar is at rest) and is dV = 
r 2 sm$drd&d<p. 

- T) is the doppler factor, which depends on the angle between 
the line of sight and the direction of motion of the wind. It is 
calculated by the expression 

£> = l — (A.2) 

T(l -BO r) 

The unit vector f is the radial unit vector in a spherical coor- 
dinate system, the center of which is located at the pulsar. O 
is a unit vector in the direction of the observer. The cosine of 
the angle between these two vectors depends on ft and tp, and 
therefore T> = D(§, tp). The observed frequency is v = Dv' . 

- We will assume that the direction to the line of sight is at po- 
lar angle f and at tp — 0. In a cartesian system of coordinates, 
then, we have: 

O = sin (x + cos (z 
whereas the radial unit vector can be expressed as: 

f — sin § cos ipx + sin § sin <py + cos §z 

which gives us: 

O ■ f — sinf sin t? cos ^ + cos £ cost? (A.3) 

- Finally, d is the distance of the object from the observer. We 
define a luminosity associated with the pulsar as: 

£ v = 4nd 2 T n F v = 4nTn j £> 2 j' W) d 3 x (A.4) 

The parameter Tn is the beam correction factor, which is 
taken to be equal to To. = (4w) _1 in the plots of the pulsar 
lightcurves and spectra, Figs. [T]|3] and |4] It is, however, taken 
to be equal to unity Ta = 1 for the L y - L sp plot in Fig. [6] 
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Appendix B: Calculation of the emission coefficient 

The emission coefficient can be calculated by multiplying the 
energy distribution of the particles by the single-electron syn- 
chrotron spectrum. This implies that the electrons are relativis- 
tic, which is, strictly speaking, not true for the whole population, 
since the distribution starts from y' — 1 . However, since the tem- 
perature of the distribution is relativistic this approximation will 
not introduce a significant error. The distribution gives the num- 
ber of electrons per unit gamma factor, per solid angle and per 
unit volume: 



b: 



while the local magnetic field in the wind frame is: 

Blc 
R 2 

D , B LC . 

B m = sin n 

v /3TR 



(B.7) 

(B.8) 
(B.9) 



dN' 



N' 



(Px'dy'dQ! 47r0# 2 (l/0) 



ties in 

/ ^T' 2 - le" r70 cosh" 2 1^1 X ' 



where » 1 is the dimensionless temperature of the distribu- 
tion in units of mc 2 as defined in the text. This distribution is 
isotropic and its density falls towards the edge of the sheet. The 
direction perpendicular to the sheet midplane is denoted by X' 
and Ax< is the width of the sheet at its local rest frame. 

Using the small argument approximation of the modified 
Bessel function Ko, Eqn.fTOl the distribution function can be ex- 
pressed as: 

dV'dy'dQ.' 20 3 ' V \a x , j 

This then has to be doubled to take account of both species in the 
sheet (electrons and positrons). The single-electron synchrotron 
spectrum is: 



V3« 



dE' 
dt'dv' 



F{x) is the well-known synchrotron function. In the above ex- 
pression, we encounter the following: 

- The field B' ± = B' sin a perpendicular to the direction that the 
photons which reach the observer travel in the wind frame, a 
is the angle between the field in the wind frame and the pho- 
tons trajectory, which can be computed through the scalar 
— > 

product B ■ O' . In order to calculate this, we need to find 
the relativistic aberration of photons which follow the line of 
sight in the lab frame. This is equivalent to calculating the 
aberration of the vector O. The field has only a radial and 
an azimuthal component, so only the aberration of these two 
components of the unit vector O are needed. These compo- 
nents are: at the point (r, §, tp) we have: 



O r = sin £ sin ft cos ip + cos ( cos § 
O v = — sin £ sin tp 



(B.2) 
(B.3) 



The corresponding primed components in the local wind 
frame are: 



O' = 



o' = 



o,-p 

1 -PO r 

O' 



r(i -po r ) 

and the angle a is 



(B.4) 
(B.5) 



a = arccos 



(6> ■ B'/B'^j 



arccos 



' O'B'+O'B' ^ 



where the superscript LC refers to quantities at the light 
cylinder. 

- The ratio X' /Ax> can be approximately converted to quanti- 
ties in the lab frame to give: 

r_ ^R-R (x,#,<p,t) 
A x , ~ 



(B.10) 



where Rq(x, ■&, ip, t) is the location of the current sheet point 
under consideration at time t. 
- The critical frequency V a is 



3eB' ± 
4mnc 



y' 2 tanh 



(B.ll) 



In the general case of the oblique rotator the temperature is 
given by the expression 



= & ± 



1 + 



\ Rq sin ■& ) 



1/3 



C(0,*) 



CSC V cot% 
, -y/l - COt 2 X COt 2 



(B.12) 



(B.13) 



where ± is the value of the temperature for the perpendic- 
ular rotator, given in Eqn.[T9l 

The function C{d-,x) is very close to unity for a large range of 
ft given a x, and deviates significantly from that value only for 
values ■& ^ n/2 ~x, or for small obliquities, as is seen in Fig. IB. II 
Therefore, the perpendicular rotator is a good approximation to 
the general problem, and this is why we have used ®± in all our 
estimates. 

The emission coefficient / has to be calculated at the re- 
tarded time 

d -f 6 

fret - t ~ - + 

c c 

where ~? is the position vector of the radiating point in the cur- 
rent sheet, d is the distance of the pulsar from the observer and t 
is the time that the observer measures. Normalizing time to l/co 
and distance to r^c the retarded time becomes 



4t - f * +R Q r-0 

rhc 



(B.14) 



(normalized time is denoted with an asterisk). The retarded time 
enters the computation via the equation for the motion of a cur- 
rent sheet: 

Ro = B (+ arccos (- cot^- cot i?) + f* et - ip) 

We replace the retarded time in the last equation from , and move 
the initial condition for time at f = d/r^c so that the new ob- 
server's time is t* d -f — d/r^c- We get then an expression for the 
radius of the current sheet as a function of observer's time, polar 
angle and azimuthal angle: 



(B.6) P (± arccos (- cot^ cot d) + t * - tp + 2&tt) 



Ro = 



i -po- 



(B.15) 
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Fig.B.l. The function C{#,x) plotted for < x < t/2,0 < § < 
n/2. Away from the sheet's edge, where # = n/2 ~x, its value is 
close to unity. 



Depending on how large or small the dot product 6 ■ r is, the 
radius Rq where the current sheet part under consideration is lo- 
cated, the radiation from which is received at time t* d , can vary 
strongly. This means that at each moment t* d the observer re- 
ceives radiation that comes from different parts of the current 
sheet that propagates in the wind, with different coordinates 
(Ro, ip). We note here that Rq is considered to be the radius 
at the midplane of the current sheet. 

Finally, the flux can be computed, according to the above, as 
follows: 



Ty 



I 



3 B c 



dvdy 1 (B.16) 



with the volume element, the function n{y'), the constant A and 
the magnetic field B' given by: 



dv = R l sin&dRd&dip 

n (Y) = y yjy 2 - \ e - y ' i& 

V3 



A = 



B' 



And 2 



2 3 



B LC (R sin# 

~RTy[—r- 



+ 1 



(B.17) 
(B.18) 

(B.19) 

(B.20) 
(B.21) 



a pulsar period, and each frequency v, produces the spectra and 
lightcurves needed. 

Finally, it should be noted that the above prescription for the 
calculation of the pulsar lightcurves is accurate as long as A < 
1/r. In the opposite case, the pulse width is not governed by F, 
but is comparable to A, something that is not taken into account 
in the calculation presented here. However, the phase-averaged 
fluxes calculated by our model are reliable even for A > 1/r. 
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The doppler factor has been given in Eq. IA.2I and is calculated 
with the help of Eq. IA.3I The strength parameter «lc and the 
density in the middle of the current sheet N' Q are given in the 
text. The expression for the flux also includes the sine of angle a 
which is calculated by Eq. IB.6l The critical frequency v' a is given 
in Eq. IB. Ill with B' ± = B' {) sin a. Finally, the above parameters 
are calculated at radius Rq with the help of t* d from Eq. IB.15l The 
purpose of the integration on R is to calculate the integral of the 
emission coefficient across the sheet, when the sheet midplane is 
at the point {Rq, §, tp). The numerical computation of the above 
integral for each point in time t* d , or equivalently each phase in 
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